Theory of Magnetic Field-Induced Bose-Einstein Condensation of Triplons in Ba 3 Cr 2 8 
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Motivated by recent experiments on BaaC^Os, a new spin-dimer compound with spin-1/2 moments of Cr 5+ 
ions, we theoretically investigate the field-induced magnetic ordering in this material in view of the Bose- 
Einstein condensation (BEC) of triplet excitations (triplons). We apply the self-consistent Hartree-Fock-Popov 
(HFP) approach to a microscopic Hamiltonian, using the realistic triplon dispersion measured in an inelastic 
neutron scattering experiment. In particular, we ask to what extent the BEC of dilute triplons near the critical 
field can explain the magnetic ordering in this material. For example, we investigate the temperature range 
where the BEC picture of triplons can be applied via the HFP approach. We also determine the temperature 
regime where a quadratic approximation of the triplon dispersion works. It is found that the strength of the 
effective repulsive interaction between triplons is much weaker in BaaC^Os than in the canonical spin-dimer 
compound TICUCI3. Small effective repulsive interaction in combination with the narrow band of triplons leads 
to higher density of triplons n CI at the critical point. It turns out that the combined effect points to a bigger 
HFP correction Un CT in BaaC^Os than in TICUCI3. Nonetheless, the HFP approach provides a reasonable 
explanation of the transverse magnetization and the specific heat data of BagC^Og. 
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I. INTRODUCTION 



Magnetic-field-induced quantum phase transitions in spin 
dimer systems have provided excellent playgrounds for the 
investigation of novel universality classes of zero tempera- 
ture quantum phase transitions. 1 These systems possess non- 
magnetic spin singlet ground states with a spin gap to triplet 
excitations (triplons). When a magnetic field H is applied, a 
quantum phase transition occurs at a critical field H c , where 
the spin gap closes and the lowest triplet excitation condenses. 
At H > H c , the average triplon density is finite and can be 
controlled by the applied magnetic field. 

The resulting ground states are determined by a delicate 
balance between the kinetic energy and the repulsive inter- 
action between triplons. 2 On one hand, if the triplon hopping 
processes are suppressed by frustration or the repulsive inter- 
action dominates, the condensed triplons may form a super- 
lattice with broken translational symmetry, leading to mag- 
netization plateaus. This is known to occur, for instance, in 
SrCu2(B03). 3 ' 4 On the other hand, if the magnetic interaction 
does not have much frustration or the kinetic energy domi- 
nates, the ground state can be described as a Bose-Einstein 
condensate (BEC) of triplons and form a homogeneous mag- 
netically ordered state. In this case, the magnetically ordered 
state at H > H c supports a staggered magnetization trans- 
verse to the field direction, leading to a canted antiferromag- 
netic state (until the system eventually becomes fully polar- 
ized as H increases). This type of behavior has been observed, 
for example, in three-dimensionally coupled spin dimer sys- 
tems TICUCI3 5 ' 6 and BaCuSi206, 7 which exhibits unconven- 
tional critical behavior. 8-11 Furthermore, recent discoveries of 
A 3 M 2 Os, 12 where A = Ba or Sr, and M = Cr or Mn, have pro- 
vided a lot of excitement for spin dimer system research, as 
these systems may represent a variety of different spin dimer 
interactions and quantum ground states. 

In this work, we present a theory of the magnetic field- 




FIG. 1: (colour online) Schematic diagram showing two neighboring 
triangular lattice planes of dimers in BaaQ^Os- Two primitive lattice 
vectors a and b are shown in the lower plane. The third primitive lat- 
tice vector c' connects central dimers in the neighboring planes. We 
set the vertical distance between neighboring planes by c/3. Here 
we use Jp to indicate the inplane nearest-neighbour interdimer cou- 
pling. Ji ( J2) denotes the nearest-neighbour (next nearest-neighbor) 
inter-plane dimer coupling. 



induced quantum phase transition discovered in BaaQ^Os, 
where Cr 5+ carries an 5=1/2 moment (3d 1 ). 13 Low temper- 
ature bulk susceptibility shows that this compound does not 
have any magnetic long-range order down to 1.5 K in the 
absense of an external magnetic field. 1415 When the exter- 
nal magentic field H reaches H c \ ~ 12 T, a field-induced 
transition to a magnetically ordered state occurs and a fully 
polarized state arises at H > H C 2 ~ 23 T. 13 In this com- 
pound, two neighboring 5=1/2 Cr 5+ ions lying along the c 
direction form a singlet dimer. In the a6-plane, these dimer 
singlets are coupled into triangular lattices, which are stacked 
along the c direction (see Fig.l). According to recent elastic 
and inelastic neutron scattering measurements, 12 BaaG^Os 
is an excellent model system for weakly coupled 5 =1/2 
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quantum spin dimers, featuring strong intradimer coupling 
of J =2.38(2) meV and weak interdimer couplings less than 
0.52(2) meV. 12 Because of the orbital degeneracy of the Cr 5+ 
ion, there is a structural transition around 70 K via a Jahn- 
Teller distortion, relieving the frustration. As a consequence, 
spatially anisotropic interdimer couplings arise. The relative 
orientations of the anisotropic interdimer couplings are de- 
scribed in Fig. 2. It was confirmed that the magnetically or- 
dered state has a commensurate and collinear transverse spin 
component for H c \ < H < H c2 . 13 This is in contrast to the 
case of Ba 3 Mn 2 8 with orbitally non-degenerate S=l Mn 5+ 
ions (3d 2 ), where the geometric frustration and single-ion 
anisotropy lead to incommensurate spiral order upon triplon 
condensation. 17 

We first consider the Heisenberg spin Hamiltonian using 
the spin exchange couplings determined by inelastic neutron 
scattering measurements. 12 Applying the bond operator for- 
malism, we obtain the dispersion of the lowest energy triplet 



excitations 
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We confirm that the Hamiltonian written in 



terms of bond operators at the quadratic level, neglecting sin- 
glet fluctuation, leads to the same triplon dispersion as deter- 
mined in experiment in Ref. 12. 

We then use the Hartree-Fock-Popov (HFP) approxima- 
tion combined with the realistic triplon dispersion to inter- 
pret two different experimental data sets, those of M. Kofu 
et al. in Ref. 13 and A. A. Aczel et al. in Ref.16. In partic- 
ular, we would like to understand to what extent Ba 3 Cr 2 08 
is a good candidate for the BEC of triplons in comparison 
to other three-dimensionally coupled spin-dimer systems such 
as TICUCI3. 20-22 Within the HFP analysis, we determine the 
effective inter-triplet repulsion U and the zero-field spin gap 
A in Ba 3 Cr20s. It turns out that the strength of the effec- 
tive repulsive interaction U between triplons in Ba 3 Cr20s is 
an order of magnitude smaller than that of TICUCI3, 21 and 
smaller than the bandwidth of the triplons as well. This sug- 
gests that the system is indeed in the regime where the ki- 
netic energy dominates. In addition, the shape of the disper- 
sion near the triplon band minimum results in the large effec- 
tive mass of triplons in Ba 3 Cr 2 08. As a result, the relation 
\H C (T) — H c (0)] oc T 3/2 in three dimensions for quadratic 
triplon dispersion works only at T < 0.1 K, while it works 
at T < 1 K in TlCuCl 3 . The HFP approach is used to de- 
scribe other experimental measurements such as the longitu- 
dinal and transverse staggered magnetizations and the heat ca- 
pacity. Despite the simplicity of the theoretical approach, the 
HFP approximation is found to explain these physical proper- 
ties even quantitatively. 

The rest of the paper is organized as follows. In Sec. II we 
use the bond operator approach to obtain the triplet dispersion 
from the microscopic Hamiltonian. Theoretical description 
of the triplon Bose-Einstein condensation within the HFP ap- 
proximation is discussed in detail in Sec. III. In Sec. IV, we 
apply the HFP approach to explain the experimental data and 
draw the phase diagram in the plane of H c versus T. The HFP 
approach is applied to describe the heat capacity in Sec.V 
and the magnetization measurements in Sec.VI. Finally, in 
Sec. VII, we summarize our result, and discuss possible limi- 
tations and extensions of the current work. 



II. TRIPLON DISPERSION VIA BOND-OPERATOR 
APPROACH 

Since the intradimer exchange interaction dominates all the 
other interdimer couplings in Ba 3 Cr20s, it is natural to take 
advantage of the bond-operator representation of the singlet 
and triplet dimer states. 23 24 This is achieved by placing one 
singlet or triplet boson on each dimer, to represent the states 
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where the quantization z-axis is taken to be the applied field 
direction. The triplet states \t m ) (m = 1, 0, —1) are chosen as 
the S z eigenstates satisfying S z \t m )=mh\t m ) . The hard-core 
constraint s^s + t^to + t\t+ + t_t- = 1 is enforced on each 
dimer, ensuring that the physical state is exactly one of the 
four above. 

The two spin operators constituting a dimer can be rewritten 
in terms of the bosonic bond operators as 



S 2a = ^ (-s f i Q - 4 s - ieapjt^t-^ 



(2) 



with a £ {x,y,z} and e the totally antisymmetric tensor. 
This form gives the correct matrix elements in the singlet- 
triplet Hilbert space. We define \t a ) triplons as eigenstates 
with S a \t a ) = so that t z = t , t x = (i_ + 1 + ) and 

t y = -j=(t- —t + ). We henceforth assume a sum over re- 
peated indices. 

Besides connecting spin operators to singlet and triplet bo- 
son operators, the bond operator formalism naturally yields 
a triplon dispersion from the microscopic Hamiltonian. We 
consider a Heisenberg spin Hamiltonian with intra-dimer cou- 
pling, and coupling between nearby dimers. We include 
nearest neighbor interactions between dimers on the same 
plane. Between adjacent planes, we include interactions 
between first and second-nearest neighboring dimers. The 
tetrahedrally-coordinated 3d 1 electron in Cr 5+ has e g orbital 
degeneracy and undergoes Jahn-Teller distortion. This gives 
rise to spatially anisotropic interdimer interactions. 12 The in- 
plane projection of the anisotropic interdimer interactions are 
described in Fig. 2. 

The intra-dimer coupling has a strength of Jo = 2.38 meV. 
Table I displays the values and directions of the couplings to 
nearby dimers depicted in detail in Fig.2. These strengths 
have been determined in Ref. 12 by fitting an RPA dispersion 
to the triplon dispersion measured by inelastic neutron scat- 
tering. 

There are four possible spin-spin interactions between two 
dimers at sites i and j. However, at quadratic order in the t 
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2a/3 + b/3 + c/3 
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= —0.15 


—a/3 + b/3 + c/3 
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Jl" 


= 0.1 


-a/3 - 2b/3 + c/3 
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= 0.04 


-4a/3 - 2b/3 + c/3 
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j" 

■J2 


= 0.1 


2a/3-2b/3 + c/3 


9 


jlll 
•J2 


= 0.09 


2a/3 + 4b/3 + c/3 



TABLE I: Interaction strength J m and relative distance Ar m for the 
nine different anisotropic interdimer couplings as described in Fig.2. 




FIG. 2: (colour online) In-plane projection of the strengths and di- 
rections of the different anisotropic interdimer couplings considered 
in this work. The couplings Jp are in-plane nearest neighbor cou- 
plings, which are represented by dashed lines. Ji (J2) is the first 
(second) nearest-neighbour interplane coupling. 



bosons, Sji • Sj-i = S j2 • S j2 = —&u • Sj2 = S i2 • Sji. 
All of the potential spin-spin Heisenberg interactions between 
two dimers thus reduce to a single effective interaction J m . 
With this, we take the Hamiltonian to be 



i i m—1 

where the intra-dimer interaction can be written as 



JqSu ' S 2 i — Jo ( :s],Si 
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We now investigate the ground state and the excitations of 
this Hamiltonian in the low temperature regime where the 
intra-dimer interaction dominates. Singlets are energetically 
favourable in this case, with macroscopic occupation, we can 
take them to be condensed. This process ignores singlet fluc- 
tuations by replacing the spin operators s\ and Sj with a com- 
plex number s, creating a uniform condensate. 

Condensing the singlets leaves a Hamiltonian describing 
the triplets. Inter-dimer interactions lead to triplon excitation 



dispersion in momentum space. Our analysis considers the 
low triplon density limit, so we retain terms up to quadratic 
order in the t operators. To impose the hardcore constraint we 
introduce a site-dependent chemical potential /Xj. This adds to 
the Hamiltonian the constraint term 
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To make a tractible analysis, the above constraint is enforced 
in a mean-field manner taking /i, = /i, on average over the 
entire lattice. In the momentum representation, this gives s 2 + 
J ( 27r )3 ^kc^ka = 1- At T = 0, the values of 11 and s are 
determined from the saddle point condition. First, / ^ \ = 



si Si 



I'iaHa 



enforces the constraint = ( 1 

as expected. Second, the condition (^p) = minimizes the 
ground-state energy with respect to the singlet density. 

Our Hamiltonian obtains a quadratic form in the momen- 
tum space: H = N^eo + T~£o + 7~L±. Here Nj, denotes the 
number of dimers on the lattice. The to triplons interact with 
themselves but not the other triplon species. They contribute 
with the quadratic Hamiltonian 



where A k = ^ - /i + £? k , with 



-Bk = ~s 2 J m cos (k ■ Ar ri 



(7) 



Ho can be diagonalized by the Bogolibov transormation 
7ko = Wk<ko + Wkf_ko> w ^ tn quasiparticle energy 25 



, - B 2 

x k D k 



H)B±. (8) 



Neither to nor 70 are subject to Zeeman splitting by the exter- 
nal field. However, the t + and t- triplons are split. Further- 
more, they interact with each other. The resultant quadratic 
Hamiltonian H± is 



W — ^2 JoSa • S j2 + JmSi,l ' Si+5r m ,l, (3) J{ ± = — 




(9) 

where h = giisH and ^ k = ^k+ *k- ^-k+ i-k-). H± is di- 
agonalized into quasiparticles 7k± with energy Wk T h, where 



7k- = «k-4- + «k-*-k+, 



(10) 



«k+ = 



Bk 



y / 2uj k (A k - UJ k ) 

A k + ^k 



--, "k+ = 



A k - ui k 



«k- = 



V^kT^lT^wk)' 
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y/2u) k (A k + UJk) 
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The 7+ triplon, with spin along the quantization axis, will be 
the focus of the HFP treatment, since it interacts with no other 
species and lowers its energy from the Zeeman splitting. Fi- 
nally, the constant part of the Hamiltonian is given by 

eo = -|Jos- 2 +M(l-5 2 )- 2^E Ak - ( 12 > 

In the limit of vanishing inter-dimer interactions (in this 
case, all J', J", J'" — > 0), the saddle-point solution gives 
s = 1 and fi = — | Jo. This limit serves as a good starting 
point in Ba 3 Cr 2 08, where J is much larger than the inter- 
dimer couplings. Furthermore, with these values of s and \i, 
the triplon dispersion matches the RPA form fitted to experi- 
mental values in Ref.12. Solving the saddle-point conditions 
with the couplings in Table I taken as bare values, we find 
s = 0.992 and /j, = — 0.775Jo, showing a high degree of 
dimerization. The bare couplings will be slightly renormal- 
ized (compared to the s = 1, fi — —3 Jo/4 case) as a re- 
sult. In particular, Jo — > Jo — A[i and J m — > s 2 J m for 
m e {1, . . . , 9}. However, we can take these renormalized 
couplings to have the values in Table I, since only the final 
dispersion will be used in our Hartree-Fock calculation. We 
furthermore assume that the dispersion is not temperature de- 
pendent within the low-temperature regime considered. 

III. HARTREE-FOCK EFFECTIVE HAMILTONIAN 

We turn our focus to the 7+ quasiparticle; being the field- 
aligned quasiparticle, it will condense with sufficiently large 
Zeeman splitting. We consider field strengths large enough 
that we may ignore the higher-energy 70 and 7_ quasiparti- 
cles. The typical splitting is of energy <?/ifci7 c (0) = 15.4 K. 
This scale is significantly larger than the highest temperature 
(around 2.7 K) where the BEC transition occurs 16 . Conse- 
quently, ignoring terms in the Hamiltonian with 70 and 7_ is 
a safe approximation to make in this external field regime. 

The Hamiltonian we take for the 6 = 7+ triplons is H = 
TLk + 7~tu, which is the sum of the kinetic and inter-triplon 
interaction terms. Here, 

Wx = ^(e fe - M )&L&k, (13) 

k 

d k,k',q 

where fi = g/J-sH — A is the chemical potential and A is 
the zero-field gap (1.37 meV from the bond-operator theory). 
£k + A is the zero-field dispersion with ek determined from 
the bond-operator theory. The quartic terms from the bond- 
operator theory give rise to interactions between the b triplons. 
Combined with the interaction from the hard-core constraint, 
this gives an approximate form for Uq. However, in the low- 
temperature limit where the excited triplons lie near the band 
minimum at Q, we may approximate rj q w Uq as a constant, 
U. The value of the interaction parameter U will be deter- 
mined from a fit to the experimental data. 



The condensate will form at the dispersion minimum Q = 
i(u + v) 13 . Note that we define the reciprocal lattice vec- 
tors u, v, and w in the conventional way. For example, u = 
2ir a ,^bxc) • W e follow the Hartree-Fock-Popov approach of 

Ref. 21 by condensing the triplons at Q: bL, 6q — > y/N^ric, 
where n c is the condensate density (the condensed boson frac- 
tion per dimer). Introducing the summation which ex- 
cludes any terms containing creation or annihilation operators 
at momentum k = Q, we decompose Hu as follows: 

nu =2N d N c+^l. j 2 + K b «] 

k,q 

+ 27T E '^L Wk'-q, (15) 

d k,k',q 

using the fact that 2Q is a reciprocal lattice vector, so that 
&-k+Q = &-(k+Q)- Performing a mean-field quadratic de- 
coupling of the quartic terms yields the following mean-field 
quadratic Hamiltonian: 

Wmf = E a + 5>~ k & k 6 k + ^ E'^-q + fe -q fe U' 

k 1 q 

(16) 

where 

E a = -im c + UN d (^--(n- n c ) 2 ^j , 

?k = £k ~ A) 

fi = gfi B H - A-2Un. (17) 

This decoupling is valid so long as the triplon densities 
6 k 6k^ are small. In the noncondensed (normal) phase, the 
Hamiltonian is already diagonalized. The triplon density is 
given by the Bose distribution function, and must be deter- 
mined self-consistently: 

/r\ 3 k 
(2^ M£ ~ k) - (18) 

In the condensed phase, we must perform another Bogoli- 
ubov transformation, which leads to the following diagonal- 
ized Hamiltonian: 

Hmf = E'^k (yt<Pk) - \ E' e ~k + Eo, (19) 

k k 

where 

Vk = u k b k + Sk&lk' ( 2 °) 

in which 

" k= /5+I' ^/SFI (2i) 
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The number of thermally excited triplons is 



d 3 fc 
(2^)3 



(22) 



For the final (p quasiparticles to be condensed at Q, they 
must be gapless. This constrains the effective chemical poten- 
tial as jx = —Un c , so that the field in the condensed phase is 
given by 



g/j, B H = A + U (2n - n c ) . 



(23) 



Between these two phases is the transition curve H C (T) 
where triplons begin to condense. The b triplons are gapped 
in the noncondensed phase. However, at the transition point 
to the condensed phase, they become gapless. This constrains 
the effective chemical potential as fl = to give the critical 
field 



H C (T) 



A 



2U 



g^B gLiB 



n a (T). 



(24) 



Since ?k = £k in this case, we determine the critical boson 
density at the transition, n cr , by the integral 



n„ = 



d 3 k 



/u(ek). 



(25) 



IV. CRITICAL DENSITY PHASE DIAGRAM 

In the HFP approach the critical field H c depends on the 
critical boson density n a linearly as shown in Eq.(24). A lin- 
ear fit of H c to n a determines the interaction parameter U 
from the slope. In addition the zero-temperature critical field 
H c (0) gives another estimate for the gap A. With a given ex- 
perimental data (H c , T), we obtain (H c ,n cr ) pairs making use 
of the Eq.(25). In the low-temperature, low-density regime 
where HFP approach is valid, we expect H c to be linear in 

™cr- 

We present three different fits of H c to n cr . These are 
based on the two different experimental data sets obtained 
from Ref.13 and Ref. 16. Lines of best fit neglect the high- 
temperature (density) regimes where H c loses linearity in n CI . 
A linear relation between H c and n a is achieved for the case 
of H parallel to the c-axis, as shown in Fig. 3 and Fig. 4. In 
Fig. 3 we describe the result of the linear fit which is ob- 
tained using the experimental data from M. Kofu et al. in 
Ref.13. From the linear relation, we estimate the interaction 
constant U = 6.5 K and zero-field spin gap A = 1.35 meV. 
The same analysis is performed using the second experimen- 
tal data set given by A. A. Aczel et al. in Ref. 16 and the 
result is shown in Fig.4. This fit gives an estimate of U s=s 8.7 
K and A = 1.34 meV, which is very close to the value ob- 
tained from the bond-operator approach. Since this experi- 
mental data shows smaller deviations from the linear fit, over 
a larger range of temperatures, we use the estimate of U w 8.7 
K and A = 1.34 meV in the following analysis. In Fig.5 we 
display the low temperature phase diagram H C (T), which is 
again obtained using the data of Ref. 16. 
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FIG. 3: The critical field H c as a function of critical density n cr 
with the applied field H parallel to the c-axis. n a (T) is calculated, 
using the HFP approach, from the given experimental temperatures. 
A linear fit of H c to n cr is performed in a low-density range. The 
experimental data is from M. Kofu et al. in Ref.13 
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FIG. 4: Same plot as in Fig.3, but the data is from A. A. Aczel et al. 
in Ref. 16. 





T(K) 

FIG. 5: Phase diagram giving the critical field H c as a function of 
temperature. Experimental data is from A. A. Aczel et al. in Ref. 16, 
with applied field H parallel to the c-axis. The solid line shows the 
theoretical result obtained from the HFP approach using the linear fit 
displayed in Fig.4. 
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FIG. 6: The critical field H c (circles) as a function of critical density 
n cr with the applied field H perpendicular to the c-axis. Here we use 
the data in Ref.13 . Inset: discrepancy between experiment and the 
best-fit line in the small n cr regime. 
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FIG. 8: Critical density n cr of BasG^Os calculated from the HFP 

3 

theory with the full dispersion (cross points). The solid line is T 2 
power-law fit coming from the simple quadratic dispersion. The 
points obtained from the full dispersion begin to deviate from the 
power-law fit around 0.06 K. 



On the other hand, the data for H perpendicular to the 
c-axis features low-temperature behaviour inconsistent with 
the general linear trend as shown in Fig. 6 and its inset. We 
think that the existence of DM interaction is one possible ex- 
planation of this low-temperature discrepancy from the lin- 
ear behavior. Further discussion on this direction is shown in 
Sec.Vn. 

Since a full-dispersion treatment successfully reproduces 
the phase diagram for the dimerized spin system TICUCI3, 21 
it is instructive to use it in comparison with the HFP approach 
for Ba 3 Cr 2 0s. Triplons in Ba 3 Cr 2 C>8 have a smaller self- 
interaction constant, of U w 8.7 K, compared to TIQ1CI3, 
which has U w 320 K. 21 However, triplon densities are signif- 
icantly higher in Ba 3 Cr 2 Os than in TICUCI3, by over an order 
of magnitude. This makes the Hartree-Fock critical field shift 
Un a greater in Ba 3 Cr 2 C>8 than in TICUCI3. In Fig. 7 we plot 
Un a in these two systems with varying temperature. Within 




0.5 1 1.5 2 2.5 3 3.5 4 
T(K) 



FIG. 7: Comparison of overall interaction energy scale Un a as a 
function of temperature between BagC^Og (solid line) and TICUCI3 
(dashed line) systems. n cr is the density of triplons at the condensate 
transition and U the inter-triplon interaction strength. 



the Hartree-Fock-Popov approach, the term —2Un acts as a 
shift in the effective chemical potential as shown in Eq.(17). 
A decrease in U will increase the effective chemical poten- 
tial, causing an increase in the triplon density as seen with 
Ba 3 Cr 2 8 . 

The shape of the dispersion affects the temperature range 

in which the power-law behavior of n cr oc is satisfied. 

At very low temperature, the quadratic approximation to the 

minimum of the dispersion becomes very accurate. As T — > 

,2 3 
0, the quadratic dispersion £k = ^ yields n cr oc T2 by 

evaluating Eq.(25) exactly 20 to give 

limn cr (T) = f (-) . (26) 

In Ba 3 Cr 2 Os, a low-temperature T ~i fit deviates from the full 
disperion critical density n cr around 0.06 K as shown in Fig. 8. 
This is an order of magnitude smaller than for TICUCI3, where 
the behaviour persists up to about 0.6 K as described 
in Fig. 9. The lower temperature scale of Ba 3 Cr 2 Os is ex- 
pected to be from the smaller triplon bandwidth, represented 
by the large effective mass near the dispersion minimum. 
From the power-law fits to Eq.(26) in Fig. 8 and Fig.9 we 
find that l/m = 1.36 K (43.6 K) for Ba 3 Cr 2 8 (TlCuCl 3 21 ) 
showing the narrower bandwidth of Ba 3 Cr 2 Os. Here we set 

h 2 /k B = 1. 

V. SPECIFIC HEAT 

We apply the HFP approach to explain the specific heat data 
measured by M. Kofu et alP To determine the magnetic con- 
tribution to the specific heat, we first find the expectation value 
of the energy per dimer. After condensing the triplons at mo- 
mentum Q, the diagonalized mean-field Hamiltonian contains 
only number operators of thermally distributed bosonic quasi- 
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FIG. 9: Same plot as Fig.8 but for the TlCuCl 3 compound. Note that 
the results obtained from the full dispersion follow the power-law fit 
up to 0.6 K. 



Non-magnetic contributions to the specific heat, such as 
phonon contribution, will not change appreciably with the ap- 
plied field. The difference Cv{H) — Cy(0) thus captures 
the heat capacity contribution from triplons. Currently, there 
exists no zero-field specific heat data, preventing proper quan- 
titative comparison. 

However, we may still make a comparison, up to an over- 
all scale difference, between the theoretical specific heat and 
experimental heat capacacity data. Fig. 10 shows the calcu- 
lated magnetic contribution with the experimentally deter- 
mined heat capacity in Ref.13. The relative scale is chosen 
to best show similarity in the peak shape for fields close to 
the zero-temperature critical field. Despite the scale differ- 
ence and non-triplon contribution, the theoretical result still 
captures the peak at the critical temperature. However, the 
drop in heat capacity is overestimated. Furthermore, it is dis- 
continuous, which can be considered as an artifact of the HFP 
approximation. 21 



particles (see Eqs. (16) and (19)). By differentiating the en- 
ergy with respect to temperature, we find the specific heat per 
dimer. In the normal phase, 
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FIG. 10: Comparison of experimental heat capacity to theoretical 
specific heat per dimer as a function of temperature. Comparisons are 
made for external fields of 12.5 T, 12.7 T, and 13 T. The experimental 
data is from Ref. 13 with the applied field H parallel to the c-axis. 



VI. MAGNETIZATION 

When H > H c , 7+ bosons condense, leading to the 
macroscopic occupation of the triplet states with the momen- 
tum corresponding to the dispersion minimum. The ground 
state wave function is then given by the coherent superposi- 
tion of the singlet and the S z =l triplet states. 1 ' 18 The density 
of the condensate determines the magnetization along the z- 
direction. In addition, the condensate supports the staggered 
magnetization which has finite transverse components (S l x ) 
and (Si V ) breaking the continuous U(l) rotation symmetry 
around the z direction. 

To determine the magnetic ordering, we begin by rewriting 
the spin operators in terms of the to, t- and t + operators. Us- 
ing the bond operator representation we obtain the following 
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relations: 



then 



(Si + S2)a = -ieal3ytj 3 ty, 

{Si-S 2 ) a = sk a +tis. 



(33) 



for a £ {x,y,z}. Due to Zeeman splitting, the t z = 70 
triplets are negligible, and we find that ((Si + £2)2:) = 
((Si + S 2 ) y ) = ((Si - S 2 ) x ) = 0. As 7 _ are similarly 
negligible, we expand the rest of the triplet t± operators in 
terms of the r y± operators. After that we ignore the terms with 
7_ because they do not contribute to expectation values. The 
average spin component per dimer along the field direction, 
which is nothing but the fraction of aligned quasiparticles n, 
is given by 



((Si + S 2 ) z ) = (tlt + -tit. 

= W, Yl (^k+7k+) («k- - «k-) = 



(34) 



Since we have condensed singlet s, the staggered component 
of the spin becomes (using u_k = Uk) 



(Silx ~ Si2x) — (t i+ + t t _ + ti + + t{ 



h.c. 



sgn(S)v / 2s 



(uq- - v Q -)$t(e iQ - ri r Q ). 



(35) 



Here Tq = \1q+/ w i tn I^qI 2 = n c . Without loss of gen- 
erality, we fix the overall phase by taking Tq to be real. Only 
the coherent condensate contributes to the transverse magneti- 
zation. Similarly, the y-component comes from the imaginary 
component of the condensate, 

(Say - S i2y ) = Sgn(J ^% Q _ - «Q_)9(e«>-«r Q ). 

(36) 



The transverse spin component thus is spatially modulated by 
the condensate wavevector Q. The transverse magnetization 
per dimer can be written as 19 



M 



v^-j^^e^'iSaz-Sa*) (37) 

(«q- - vq-) e<Q ri cos (Q • r *)r Q 

i 

V2N d to--«Q-) r * (38) 
The square of the transverse magnetization per Cr 5+ ion is 



N d 

sgn(B Q )V 2s 
sgn(B Q )s 



Ml = (gn B M : 



•HI) 
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2lj, 
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8N d 2lo q (A q +uj q ) 

Q _2 Jo9 2 2 „ m 

= s Uc ~8A^b- (39) 



Having neglected the 70 and 7_ triplons, we estimate s 2 = 
1 — n, using the overall triplet boson constraint. The total 
and condensed triplet densities, n and n c , are determined by 
solving Eq.(22) and Eq.(23) self-consistently. 

The transverse magnetization has been measured by the 
elastic neutron scattering experiments. 13 The applied field 
is perpendicular to the c-axis. Fig. 11 compares theoretical 
squared perpendicular magnetization at T = 0.2 K to the 
experimental results. 13 Deviation from the experiment occurs 




13.5 



FIG. 1 1 : Comparison of experimental and theoretical perpendicular 
magnetization squared at T = 0.2 K. Perpendicular magnetization 
is defined in Eq.(37). The experimental data is from Ref.13, with the 
applied field H perpendicular to the c-axis. 

most prominently in the critical field. This is caused by dis- 
crepancy between the linear fit H c oc n cr (T) and the experi- 
mental critical field H C (T). However, the shape of the mag- 
netization curve past the critical field is properly reproduced. 
This can be seen in Fig. 12, where the theoretical result has 
been translated to match the experimental critical field. The 
resulting shape matches over the entire range of fields, with 
the theoretical magnetization larger by a factor of 1.13. The 
magnetization also jumps slightly at the critical field. Like 
the discontinuity in specific heat, this is an artifact of the HFP 
treatment. 29 

The parallel magnetization has been measured as a func- 
tion of applied field (both parallel and perpendicular to the 
c-axis) at the condensate transition. 13 Fig. 13 gives the HFP 
result for H perpendicular to the c-axis, with the magnetiza- 
tion per dimer My = |^sn. The saturation field, where all 
spins are aligned with the field, is severely underestimated by 
the HFP result of 18 T. Experimentally it is found between 
23 T and 24 T, from the derivative of magnetization |^. 13 
This happens because the triplon density in the HFP approach 
grows too quickly with increasing field. 
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FIG. 12: Comparison of the shape of theoretical and experimental 
perpendicular magnetization curves as in Fig. 1 1 . Theoretical result 
has been translated by 0.07 T to match the experimental critical-field 
behaviour. The scale of the theoretical result is 1.13 larger than that 
of the experimental data. 




FIG. 13: Parallel magnetization per Cr atom, as a function of applied 
field perpendicular to the c-axis, at T = 0.4 K. The result is found 
from the triplon density as determined by the HFP approach. Critical 
field H c and saturation field H s are indicated. 



VII. DISCUSSION 

We have applied the HFP approach to understand the triplon 
BEC in Ba 3 Cr 2 08, using the full dispersion of the triplons 
measured in the recent neutron scattering experiments (which 
is recast in the form of a bond-operator representation of a 
Heisenberg model). We investigated the temperature range 
where the HFP approach is valid with the full dispersion, and 
also locates the temperature where the quadratic approxima- 
tion of the dispersion breaks down. Using this approach, we 
computed the transverse magnetization and specific heat that 
are favorably compared to available experimental data. Our 
results show that the BEC picture overall works reasonably 
well for Ba 3 Cr 2 08. 

In the much-studied three-dimensionally-coupled spin- 
dimer system TlCuCl 3 , the triplon band width W ~ 87 K 
and the effective interaction U <~ 340 K within the HFP 



analysis. 21 In contrast, our analysis leads to W ~ 21 K and 
U ~ 8.7 K in Ba 3 Cr 2 8 . Thus it may appear that the HFP 
would work better for Ba 3 Cr 2 Os because of smaller U /W. 
On the other hand, smaller U in Ba 3 Cr 2 Os results in a larger 
critical triplon density n cr <~ 0.1 compared to n cr ~ 0.002 
in TlCuCl 3 , making the dilute triplon density approximation 
less valid. In the end, the combined effect in the form of the 
HFP correction to the critical field H C (T), Un cr , turns out 
to be bigger for the case of Ba 3 Cr 2 Os- This means that the 
temperature range where the HFP approach is valid is more 
limited in the case of Ba 3 Cr 2 Os- Indeed, it is found that the 
HFP works up to 8K in TlCuCl 3 while it fits the data up to 
2K at best in Ba 3 Cr 2 O s . 

The triplon dispersion in Ba 3 Cr 2 Os is flatter (or the ef- 
fective mass is larger) compared to TlCuCl 3 , which leads 
to a smaller window of temperatures where the quadratic- 
dispersion approximation is valid. This is seen in how the 
relation [H C (T) - H c (0)] oc T 3/2 reproduces the phase di- 
agram for T < 0.1 K for Ba 3 Cr 2 O s , but for T < 1 K in 
TlCuCl 3 . 

A useful way to improve the HFP results may be to in- 
troduce the hard-core constraint among the triplons. The so- 
called Bruckner bond operator approach 26 27 achieves this by 
introducing an infinite on-site triplon repulsion by 



He = V^tjjlpUaUp 



(40) 



as V — > oo. In the low-density limit, this hard-core interaction 
may be treated exactly by a summation of ladder diagrams at 
the one-loop level in the self-energy. This approach, when 
generalized to finite temperature, should lower the triplon den- 
sities and extend the region where low-density approxima- 
tions are valid. This may be a useful future extension of our 
work. 

Recent ESR measurements indicate the existence of siglet- 
triplet mixing in the ground state of Ba 3 Cr 2 Os. 13 In the 
ground state, singlets mix with to for H _L c, and with 
t± for H || c. This mixing points to the existence of a 
Dzyaloshinsky-Moriya (DM) interaction of the form • S j x 
Sj, with Djj perpendicular to the c-axis. Since it breaks the 
U(l) symmetry of the Heisenberg Hamiltonian, the system is 
no longer described by a BEC transition. The result is that 
triplons are gapped and always condensed to some extent, 
turning the transition into a crossover region. 28 ' 29 Below the 
temperature scale of the DM interaction, then, we expect that 
a simple BEC picture of triplons will no longer be sufficient. 
This could explain, for instance, the nonlinearity of critical 
field H c in critical density n cr at low temperatures for H per- 
pendicular to the c-axis. An understanding of the magnitude 
and direction of the vector is important for a proper and 
full description of Ba 3 Cr 2 Os, especially at very low tempera- 
tures. 
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